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When an electron interacts with phonons, the electron can exhibit either free electron-like or 
polaron-like properties. The latter tends to occur for very strong coupling, and results in a phonon 
cloud accompanying the electron as it moves, thus raising its mass considerably. We summarize 
this behaviour for the Holstein model in one, two and three dimensions, and note that the crossover 
occurs for fairly low coupling strengths compared to those attributed to real materials exhibiting 
conventional superconductivity. 



I. INTRODUCTION 

The conventional theory of superconductivity consists 
of two distinct pieces. The first describes the formal- 
ism of pairing, premised on the idea that two electrons 
effectively attract one another; this pairing is well de- 
scribed by the Bardccn-Coopcr-Schricffcr (BCS) wave 
function, 1 and continues to be used extensively to de- 
scribe almost all known superconductors. The second 
piece concerns the origin of the attractive interaction, 
suggested originally by Frohlich 2 - and Bardeen^ to be the 
electron-phonon interaction. In the BCS calculation the 
additional piece concerning the origin (i.e. mechanism) 
of pairing enters only in that a simple (though confusing) 
cutoff of order the Debye (phonon) frequency is included 
in the attractive interaction. It is confusing because it 
is added in momentum space, whereas the physical cut- 
off is in frequency space, to reflect the retardation ef- 
fects expected from relatively sluggish phonons interact- 
ing with nimble electrons. The cutoff occurs naturally in 
frequency space in the later formulation of Eliashberg 4 
and others 

A very solid and coherent case has been made for this 
description of phonon-mediated mechanism of supercon- 
ductivity for simple metals like Lead.— However, even 
though Eliashberg theory is often referred to as a mi- 
croscopic theory of superconductivity (and, indeed, it 
is more microscopic than BCS theory), it retains some 
elements of a phenomenology. This occurs in two dis- 
tinct areas in the theory; the first concerns the treat- 
ment of direct Coulomb interactions, which in both BCS 
and Eliashberg treatments tends to be modelled through 
the use of a single (renormalizcd) parameter, fi* , and the 
second concerns the treatment of the phonons. Tradi- 
tionally, the phonons are taken from experiment, notably 
tunneling^ or neutron scattering^ In either case all the 
"calculations" concerning the phonons have been done 
by nature, and handed over to the theorist to incorpo- 
rate into Eliashberg theory. Notable exceptions are, for 
example, treatments in Ref. (fl2r ) and more recently in 
Ref. (jl3T ). where an attempt is made to compute modifi- 
cation to the properties of the phonons alongside those of 
the electrons, due to strong electron-phonon interactions. 

As an example we show in Fig. 1 a plot of the 
pairing susceptibility (a measure of the superconduct- 



ing critical temperature) vs. electron filling, using a 
Migdal-Eliashberg calculation with and without phonon 
renormalization. 14 When phonon renormalization is in- 
cluded, the susceptibility (here plotted for an inverse 
temperature, /3 = 6), decreases, particularly as one ap- 
proaches half filling. In fact our calculations stop because 
a charge density wave instability occurs for a density less 
than half filling. This makes apparent that phonon renor- 
malization has an important but deleterious effect on su- 
perconductivity, at least for the model considered in Ref. 
(|14l ). This is an important point which will have to be 
reconsidered in future work. For the present, however, 
we wish to focus on a simpler model, that of a single 
electron interacting with many phonons. This object is 
usually referred to as a polaron (though not in the lit- 
erature pertaining to superconductivity^), and we re- 
gard the polaron as the basic building block, from which 
a theory of interacting electrons and phonons must be 
built. This is not how the theory of superconductivity 
progressed historically. A weak coupling approach was 
assumed at the onset; this means that a (large) Fermi 
sea is taken as 'given', and then Migdal-like 19 arguments 
are adopted to essentially use weak coupling perturbation 
theory3 

In this work we wish to illustrate the polaronic na- 
ture of the electron when coupled to phonons, using 
the Holstein models We will summarize results in one, 
two, and three dimensions, adopting various strategies: 
weak and strong coupling perturbation theor y 22 i 23 the 
adiabatic approximation^™ and an exact diagonalization 
method which incorporates Bloch's Theorem which we 
refer to as the Trugman methodi— ~— Very little work 
has been done beyond the Holstein model. When a non- 
local interaction is used, the coupling in momentum space 
acquires momentum dependence, making all techniques 
more difficult. Very recently we have succeeded in ob- 
taining analytical results for the polaron in the weak 
coupling limit 28 of the Barisic, Labbe, and Friedel22- 
Su-Schrieffer-Heeger— i (BLF-SSH) model, which models 
acoustic phonons coupled to electrons through their itin- 
eracy. Other, more complicated models have also been 
proposed. For example in Ref. ([31T ) a longer range elec- 
tron phonon interaction is proposed, and Jahn- Teller- 
type distortions in the presence of repulsive Coulomb in- 
teractions also lead to polaron formation^ Very recently, 
Innocenti et a/i— ■ have proposed a two electron band see- 
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FIG. 1. (Color online) Pairing susceptibility vs. electron fill- 
ing for a fixed bare coupling strength, Ao = 2, and a fixed 
temperature, given by /3 = l/(ksT) = 6. Results for two dif- 
ferent frequencies are shown; usual Migdal calculations of the 
electronic properties with no changes to the phonons are given 
by the curves, while results that allow the phonons to renor- 
malize due to the electron phonon coupling are given by the 
symbols. Clearly allowing for phonon renormalization leads 
to a suppression of the pairing susceptibility, and therefore 
superconductivity. 



nario, where tuning of the chemical potential can lead 
to polaron formation; a companion proposal 3 — also sug- 
gests that so-called 'dopant self-organization' can lead to 
trapping of polarons, resulting in some kind of inhomo- 
geneous phase. These are all interesting possibilities, but 
here we focus on rigorous results for the simplest model 
describing polarons. 

In the next section we describe the various techniques 
used to attack the Holstein polaron, followed by results 
for two specific properties — the ground state energy 
and the effective mass. It is seen, for example, that the 
electron is always polaronic in one dimension. More pre- 
cisely, this statement means the following: as the phonon 
frequency is taken to be smaller and smaller (this is neces- 
sarily a process taken as a limit) more and more phonons 
are required to properly describe the ground state, irre- 
spective of the coupling strength. In the past researchers 
have referred to this as a self-localized (or self-trapped) 
state; this is a misnomer, as no such violation of Bloch's 
Theorem occurs. Partly this nomenclature has arisen 
from variational studies (and the adiabatic limit, where 
the phonon frequency lue is taken to be zero at the on- 
set), where then a localized state does arise. The result 
is that the entity of interest emerging out of the ground 
state is hardly recognizable as an electron; rather it con- 
tains a considerable admixture of phonons, and hence is 
renamed a polaron. After presenting results in two and 
three dimensions, where a regime of weak coupling exists 
in which the electron definitely retains free-electron-like ( 
and therefore 'non-polaron-like') character for weak cou- 



pling strengths, we will conclude with a summary. 

II. METHODS 

Perhaps the simplest starting point for all studies such 
as this one is the weak coupling approximation. This is 
well-defined with straightforward pcrturbativc methods, 
or one can adopt the language of Green functions, where 
the electron propagator, G(k,uj), can be written in mo- 
mentum (fc) and frequency (w) space in terms of a self 
energy, E(fc,w): 

Z(k,w + iS) = -^2\g(k,k')\ 2 G (k',oj + iS-uj(k-k')), 

k' 

(1) 

where Go(k,oj + iS) = \w + iS — e^l is the non- 
interacting electron retarded propagator, e& (w q ) gives 
the dispersion relation for the non-interacting electrons 
(phonons), respectively, and g(k,k') describes any possi- 
ble momentum dependence in the electron-phonon cou- 
pling. Note that the infinitesimal imaginary part to the 
frequency, i5, is necessary for the analyticity of the Green 
function; here we use the retarded functions (positive 
imaginary part). This is written quite generally, and for 
any dimension. Here, we wish to discuss specifically the 
Holstein model, with Hamiltonian written in momentum 
space as 

H = ^ tkc\c h + u E ^2 a\a q 

k q 
VJV k,k> 

where uje is the Einstein phonon frequency and g is 
the (dimensionless) electron-phonon coupling, and in this 
case is clearly a constant. We will tend to use the (also 
dimensionless) parameter A, which will have a different 
definition depending on dimension. In one dimension, 
we use A = g 2 w E /(W/2) = g 2 w E /(2t), where W = At is 
the electron bandwidth written in terms of the nearest 
neighbour hopping parameter, t. Hereafter t will be set 
to unity. 

In one and two dimensions, the momentum sums can 
be evaluated analytically. The results are 

H 1B (w + tS)= _ , (3) 

so that the ground state energy is 

and the effective mass m* is given by 

X n~ 1 4- ^m. 
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FIG. 2. The ground state energy and effective mass vs electron phonon coupling strength for the Holstein model. The energy 
((a) (c) (e)) and effective mass ((b) (d) (f)) are plotted in one, two, and three dimensions, respectively. Note the very abrupt 
crossover to polaron-like behaviour in 2D and 3D, both at relatively low coupling strengths. The strong coupling limit agrees 
very well with the adiabatic limit in all three dimensions. Also note that in (a), (c) and (e) there are three curves that are 
indistinguishable for A > 1.5, showing that the strong coupling and adiabatic regime is readily achieved for moderate A and 
non-zero loe- 
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In two dimensions we obtain^ 



S 2 d(w + id) = K 

2 u — loe 



4f 



10 — uje 



(6) 



the strong coupling limit, converging results there, and 
using these results as subsequent 'seeds' as we span the 
coupling by lowering the coupling strength a little at a 
time. 



where K(x) = J^ 2 d9 



is the complete Elliptic 



\J l — x sin 2 8 

integral of the first kind. This leads to a ground state 
energy, which, in weak coupling, is: 



E Q = -At 1 + - 



Awe 



-K 



4 t l + w B /(4t) L(l + w B /(4t)) _ 
We can take the derivative of Eq. © to obtain: 



(7) 



m*/m = 1 + 



A 



2 l + w B /(8t) 



E 



(l+ W£ /(4i)) 2 



(8) 



where E(x) = J^ 2 dO\J 1 — x sin 2 8 is the complete El- 
liptic integral of the second kind. Three dimensions is 
straightforward, but requires numerical integration. 

The calculation in the strong coupling limit is outlined 
in Ref. (jH). The unperturbed state is: 



IV>> 



|0>, 



(9) 



where the sum is over all lattice sites^ Straightforward 
application of perturbation theory yields a dispersion, al- 
beit with exponentially suppressed bandwidth. This con- 
stitutes an estimate for the effective mass, which is far 
too high, and will be omitted from the plots. However, 
second order perturbation theory gives a substantive cor- 
rection to the ground state energies; these are given, for 
large coupling values, as 



Eq « -UJE9 - 



dt 2 

uj E g 2 



(10) 



where d is the dimension of the space. In the parts of 
Fig. 2 showing the ground state energy, the curve display- 
ing Eq. (|10[) is indistinguishable from the adiabatic limit 
for A Z 1. Adiabatic limit results are determined through 
a decoupling procedure, whereby the ion displacements 
are treated as c-numbersi 22 ' 24 

Finally, exact results are obtained using the Trugman 
method ) 26 ' 27 We are able to obtain converged results for 
very small values of the dimensionless parameter uiE/t by 
using Eq. ^ as a 'seed' state for the state construction. 
Iteration is considerably improved as well by starting in 



III. RESULTS AND DISCUSSION 

Results are shown in Fig. 2 for the one dimensional (a) 
and (b), two dimensional (c) and (d), and three dimen- 
sional (e) and (f ) ground state energy and effective mass, 
respectively, vs. electron phonon coupling strength, for a 
variety of phonon frequencies. Note that we use wg/i and 
A as the two variables to characterize the adiabaticity and 
coupling strength, respectively, of the electron phonon 
system. Here A = ujEg 2 /(W/2) in one and three dimen- 
sions, where W is the bandwidth (here, for a tight binding 
model with nearest neighbour hopping only, W — At in 
ID and W = 12t in 3D. In two dimensions we defined A 
slightly differently^ as A = uJEg 2 /(27ri); this is because 
the non-interacting electron density of states (EDOS) is 
a constant at the bottom of the band, and so we adopted 
this value to multiply the bare value, Ao, rather than the 
average EDOS, which is what we did in ID and 3D. With 
this definition, for example, the effective mass becomes 
1 + A/2 in 2D in the weak coupling adiabatic limits 

The main results can be summarized as follows: in 
ID the electron is always polaron-like. Note in partic- 
ular that as the phonon frequency decreases the effec- 
tive mass increases dramatically (Fig. 2(b)), and essen- 
tially diverges in the adiabatic limit, even for infinitesimal 
coupling strength. The results for the ground state go 
smoothly over to the adiabatic limit, as is apparent from 
the disappearance of two of the finite frequency curves 
into the adiabatic curve in Fig. 2(a). In 2D and 3D 
there are clear delineations of free electron-like behaviour 
at weak coupling and polaron-like behaviour at interme- 
diate and strong coupling. However all the results are 
'smooth' for a non-zero phonon frequency — only when 
an adiabatic calculation is performed (with the phonon 
frequency set equal to zero at the onset) does the curve 
display a kink (and therefore a transition). So no phase 
transition occurs as long as =/= 0. The remarkable 
result is that the crossover to polaron-like behaviour oc- 
curs at such an intermediate coupling strength, 36 well 
below the value of coupling strength normally assigned 
to conventional superconductors (even in 2D, if the av- 
erage EDOS is used in the definition of A, the crossover 
coupling strength only moves up to about 0.86, just a 
little lower than the value in 3D). Ongoing work is inves- 
tigating this crossover behaviour when acoustic modes 
are involved^ 
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